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OS ■ Complex formalism of Riemann - Silberstein - Majorana - Oppenheimer in Maxwell 

electrodynamics is extended to the case of arbitrary pseudo-Riemannian space - time 
| in accordance with the tetrad recipe of Tetrode - Weyl - Fock - Ivanenko. In this 

approach, the Maxwell equations are solved exactly on the background of static 
cosmological Einstein model, parameterized by special cylindrical coordinates and 
realized as a Riemann space of constant positive curvature. A discrete frequency 
spectrum for electromagnetic modes depending on the curvature radius of space and 
three parameters is found, and corresponding basis electromagnetic solutions have 
O ,1 been constructed explicitly. In the case of elliptical model a part of the constructed 

solutions should be rejected by continuity considerations. 

Similar treatment is given for Maxwell equations in hyperbolic Lobachevsky 
model, the complete basis of electromagnetic solutions in corresponding cylindrical 
coordinates has been constructed as well, no quantization of frequencies of electro- 
magnetic modes arises. 

> , 

^ ■ 1 Introduction 

00 

It is well-known that Special Relativity arose from investigation of symmetry properties of the 
Maxwell equations with respect to inertial motion of the reference frame: Lorentz [1], Poincare 
[2], Einstein [3]. Naturally it was electromagnetic field that was the first and principal object 
for the Special Relativity: Minkowski [4], Silberstein [5], Marcolongo [7], Bateman [8]. In 
1931 Majorana [10] and Oppenheimer [9] proposed to consider classical Maxwell equations as 
a quantum photon equations. In this context they introduced 3-vector function obeying Dirac- 
like massless wave equation. It turned out that much earlier in 1907 the same mathematical 
translation of classical Maxwell theory was performed by Silberstein [5]; besides, he noted himself 
that the same approach was used earlier by Riemann [6]. That history was much forgotten, and 
many years this complex approach to electrodynamics was connected mainly with Majorana and 
Oppenheimer. Historicaljustice was rendered by Bialynicki-Birula [11], see also in [12-17]. 

In the present paper Qj we use the complex formalism of Riemann - Silberstein - Majorana - 
Oppenheimer in Maxwell electrodynamics extended to the case of arbitrary pseudo-Riemannian 
space - time in accordance with the tetrad recipe of Tetrode - Weyl - Fock - Ivanenko (for 
more detail, see |19j). In this approach the Maxwell equations are solved exactly on the back- 
ground of simplest static cosmological models, Riemann and Lobachevsky spaces of constant 



X 



1 It is an extended version of the paper: Bogush A. A., Krylov G.G., Ovsiyuk E.M., Red'kov V.M., Maxwell elec- 
trodynamics in complex form, solutions with cylindric symmetry in Riemann space of constant positive curvature. 
Doklady of the National Academy of Sciences of Belarus. 2009 (in press) . 
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curvature parameterized by cylindric coordinate (many years ago these coordinates were used 
by Schrodinger in his book [19]; systematic treatment of coordinate systems in Riemann and 
Lobachevsky spaces was given by Olevsky [21]). In the case of compact Riemann model a dis- 
crete frequency spectrum for electromagnetic modes depending on the curvature radius of space 
is found. In the case of hyperbolic Lobachevsky model no discrete spectrum for frequencies of 
electromagnetic modes arises. 



2 Complex matrix form of Maxwell equations 

Let us start with the Maxwell equations in vacuum: 



div cB 
divE = L 



0, 



rot E 



dcB 



dct ' 
<9E 

rot cB = uocJ + — — . 

act 



(2.1) 



With notation j a = (p, J/c) , c 2 = 1/eo/^o they read 



div cB = 



P 



rot E 



dcB 



div E = — 

eo 

or in the explicit component form 

dicB 1 + d 2 cB 2 + d 3 cB 3 = 
dsE 1 - d 1 E 3 + d cB 2 = 
d x E l + 8 2 E 2 + 8 3 E 3 = f/e , 
d 3 cB l - d lC B 3 - d E 2 = j 2 /e , 



dct 

^ j <9E 

rot cB = h — — 

eo act 



(2.2) 



0, 
0, 



8 2 E 3 - 8 3 E 2 + OocB 1 = , 
diE 2 - d 2 E x + d cB 3 = , 
8 2 cB 3 -8 3 cB 2 -OqE 1 =j 1 /e , 
dicB 2 - d 2 cB x - d E 3 = j 3 /e . 



(2.3) 



With the use of complex 3-vector field ip k = E k + icB k eqs. (|2.3p can be combined into 

d^ 1 + d 2 v° + a 3 ^ 3 = j°/e , 

-idoip 1 + (d 2 ^ 3 - a 3 V 2 ) = if/eo , 

-id^ 2 + (d^ 1 - d^ 3 ) = i f/e , 
-id ^ 3 + (d^ 2 - d^ 1 ) = i j 3 /e . 

These four equations can be presented in the matrix form: 



(2.4) 



-id 



a 





1 



a 2 1 
a 3 1 



+ di 



3 



bo 


1 










h 























-1 


h 





1 







do 










1 


di 







-1 





d 2 


1 










<h 














+ d 2 



co 





1 





ci 








1 


C2 











C3 


-1 









+ 





V> 3 



l 

eo 



i f 
i f 
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There arise four matrices (including arbitrary numerical parameters) 



a 







ia°d + oPdj)V = 


J, 






1l> 

tp 
















1p 


bo 


1 







CO 





1 





bi 










Cl 








1 


b 2 





-1 




C2 











b 3 





1 




C3 


-1 









a 



a 



d 



d 2 



d 3 



a 

oi 1 

a 2 1 

a 3 1 

1 

-1 







and 



0\2 



(« U ) 



aoao 

a\ao + a\ 10 

a 2 a + a 2 1 

a 3 a + a 3 1 



Let us require 

(a ) 2 = +/, 
the most simple solution is 

a = ±1, a. 

In the same manner 

(a 1 ) 2 = 

we get 



a a = 1 , aia + a\ , a 2 a + a 2 , a 3 a + a 3 ; 



0, 



a 



±10 
10 
10 
1 



(« u ) 



0\2 



+1 



b 2 + 6i 6 

&2&0 - &3 &2 
- ^2 &3 







-1 





(a 1 )* = -/ , 



6 = , 6i = -1 , & 2 = , 6 3 = , 



a 














-1 




















-1 








1 






Analogously 



2\2 



(<0 



CO Co + c 2 





CO 





ClC + c 3 


-1 


Cl 





C2C0 





C2 





c 3 c - Cl 





C3 


-1 



(2.5) 



(2.6) 



(2.7) 
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that is 



c = , ci = , c 2 



And finally 



3\2 



that is 



d = , di = , d 2 = , d 3 
Consider their products: 



-aV 







n 
u 


n 
u 


1 
1 


n 
u 


= o, 


a 2 = 



— i 




n 
u 




n 
u 


1 
u 









-1 








dod + ^3 








do 






dido - d 2 
d 2 d + (ii 


-1 




-1 




d 2 




-/ , 


d 3 d 








d 3 










u 


L) 


U 


1 


= -1, 


a 3 = 







1 


-1 












-1 











aV = 


-a 3 a 2 


= a 


l 

1 


Q 


V 



2\2 



(<0 



-I . 



(a 3 ) 2 = -/ . (2.9) 



a 



Consider the product a a u . 



±1, aV 






A; 








-1 




















-1 








1 






1 

-At 



1 



If k = +1, we will have the most simple commutation rule: 

a° = I, 

Thus, the Maxwell matrix equation looks 



do? 



a 



(-id + cPdj)V = J 





V> 3 



J = 



a 



1 

-10 



1 



a 1 a 2 



-a 2 a l 



a 









1 














1 


-1 














-1 








(a 1 


) 2 = 


-I 


5 



a 



a 



a 2 a 3 



-a 3 a 2 



[a 
a 1 



1\2 



1 ij 1 

e * f 
ij 3 

1 
0-10 
1 
-10 

(a 1 ) 2 = " 



1 3 

—a a = a 



(2.8) 



(2.10) 



(2.11) 



(2.12) 
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3 Maxwell matrix equation in Riemannian space 

Maxwell equation 



(a do + a j dj) = J 





E + icB 



J 



eo 



-il 
P 



(3.1) 



can be extended to an arbitrary Riemannian space-time in accordance with general tetrad recipe 
of Tetrode- Weyl-Fock-Ivanenko (see |18j): 



a p (x) [ d p + A p {x) } ¥(s) = J(x) , 



af(x) = a c e[ c) {x) , A p (x) = p ab ef a) V p e (n)/9 



(3.2) 



where (x) is a tetrad; j ab stands for generators for complex vector representation of orthogonal 
group SO(3.C); V p denotes a covariant derivative. Eq. (|3.ip can be rewritten differently 



1 



» c (e1 c) d p + -j ab labc )^ = J(x), 



(3.3) 



'(c) 

with the use of Ricci rotation coefncients7;, ac = — 7 a f, c = — e (b)/3;a e f a -\ e ? c ) ■ 

Eq. (|3.ip is invariant under gauge transformations of the local Lorentz group (see |18j ) 

\P'(x) = S(z)¥(s) , S(x) e SO(3.C) loc , 

e '(a)a( X ) = L a( X ) e {b)a {x) , 

aP(x) [d p + A p (x)} *(x) = J(x) , 
a p {x) [d p + A p (x)] #'(x) = J'{x) . (3.4) 



4 Tetrad explicit form of Maxwell matrix equation 



Matrix Maxwell equation (|3.3p can be written as 



" i ( e ( ) d P + \j ab laW + a k ( e[ k) dp + \j ab labk )* = ■/(*) . 



Taking into account the identities 



1 



-f^abO = [51(7230 + «7010) + 52(7310 + *702o) + S3(7l20 + *7030 ) ] , 
^j ab 7abk = [Sl(723fc + ilOlk) + S 2 (l31k + «702fc) + S 3 (ji 2 k + «703fc) ] 



and using notation 



(0) 9 P = a (0) > e '(fc) u P = . 



c\ = dt, 



a = 0,1,2,3 , 



(701a,702a,703a) = V a , (723a, 731a, 7l2a) = Pa 



(4.1) 



(4.2) 



(4.3) 
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eq. (|4.1|) is reduced to 

- i [ 9 (0) + s( Po + iv ) ] * + a fc [ d (fc) + s( Pfe + iv fe ) ] * = J(x) , 

or 

( a k + sv + a fc sp fe ) 
-i ( a (0) + sp - a fc sv fc ) 
where s, stands for the generators 





E + icB 






1 


P 


E + icB 


eo 


i 3 






































































1 










-1 





si = 











-1 


, s 2 = 














, S3 = 





1 
















1 










-1 























(4.4) 



5 Cylindric coordinates and tetrad in spherical space S3 

Let us consider the Maxwell equation in the cylindric coordinates and tetrad in spherical space 
S3: 



n\ = sin r cos <p , n 2 = sin r sin <p , 713 = cos r sin z , = cos r cos z 
dS 2 = dt 2 — dr 2 — sin 2 r d(j) 2 — cos 2 r dz 2 , x a = (t, r, (f>, z) 



10 

10 

sm -1 r 




cos 1 r 



100 

0-10 

-sinr 

— cos r 



where (r, (f), z) run within 

P G [0, +7T/2] , (f) G [-7T, +7T] , 2 G [-7T, +7r] . 

Christoffel symbols are given by r° CT = , r oo = , r j = and 



jk 





— sin r cos r 
sin r cos r 



1 jfc 



COST 

smr 

COST q q 

sin r 



jk 





For covariant derivatives of tetrad vectors we get 
A 








sin r 
cos r 









dA ? V A 



dei 



e (0)/3;° - -g^T ~ 1 e (0)0 ~ U > 



(5.1) 



(5.2) 



6 



3(l)/3;a 



a/3 



de 



2(2)/3;a 



(2)/3 



dx c 



a/3 e (2)0 



3(3)#« 



(3)0 



r a/3 e (3)z 


























— sin r cos r 










E 




















cos r 












































- 


sinr 



























It remains to find Ricci rotation coefficients: 

7aw = e (a ^ e {b) p. t ej 0) = , 

P 



labl 
7a63 



2(a) e (b)a;r 6(1) , 

z 

2 (a ) e (6 ) / g ;z e (3) ; 



lab2 - e (a) e {m<j> e (2) , 
from which it follow 

7011 = 7021 = 7031 = , 7012 = 7022 = 7032 = , 7oi3 = 7023 = 7033 = 



7231 = , 7311 = , 

7232 = , 7312 = , 7122 

sinr 

7233 = , 7313 

Taking into account the identities 



7121 = , 
cos r 



cos r 



70) w p 



8, 



1 



(2) 



V 

v 2 = 

Po 

P2 = (7232 , 7312 , 7122) = 



(2) P v , gmr 
= (7010 , 7020 , 7030 ) = 

(70120,7022,7032) = 

= (7230,7310,7120) = 
cos r N 

(0,0, 



(0) 
rh 



sinr 

7123 = . 

dt , 



71) dp — ®t 



(i) 



,8„ = d, 



1 



(3) 



Or 



-0, 



sinr 



P3 



(3) P \ cosr 

v i = (7011, 7021, 7031 ) = 

V3 = (7013,7023,7033) = 

Pi = (7231,7311,7121) = 
sinr 

(7233,7313,7123) = (0, ,0) 

cosr 



in the absence of an external source eq. (|16,3j) reads 

( -id t + a l d r + c? — c^ + a 3 J- dz + a 2 S 3 — + a 3S 2 ^) 
\ sinr cosr smr cosr / 



6 Separation of variables in S3, solutions at m = 

Wave Maxwell operator from (|5.4p commutes with the following ones: idt, id^, id z . Therefore, 
for a field function we get a substitution 





E + icB 



Correspondingly, eq. (|5.4p reads 





/i(r) 
/2(r) 



(6.1) 



-to + a 



d 



im 



dr sin r 



a; + 



cos r 



a 3 + 



cosr 2 smr 3 
— a D3 H a 02 



smr 



cos r 





AW 

/ 2 (r) 
/s(r) 



After simple calculation we get the following radial system: 
t d cos r sin r im ik 

+ ~. )./l + ~. J2 + — J3 



dr sin r cos r 



-uf 2 



sin r cos r 

-w/i f 2 + -. — J3 

cos r sin r 

/ d sinr^ ifc 

J3 H h 







v <ir cos r y 



f , / d cosr 



cos r 
im 



0. 

/i = 0. 



(6.2) 



" (ir sin r ' sin r 

Let us consider first the case m = 0, when the radial equations become more simple 



(6.3) 



d 



+ 



cos r sin r 



dr sin r cos r 



)fl + 



ik 



cos r 



/ 3 = 



ik 



11 



-ujf 2 



uj cosr 
sinr, „ ik „ 

^ h + — /l = 

dr cos r cos r 



/3 



1 ( d 



cos r 



<ir sin r 

Using 2nd and 4th equations, from the first one it follows 

sinr . 



2 ■ 



(6.4) 



d cos r 
+ 



ik ik 1 d cosr 

/2 + —rrz — \~r + ~ — ) h 







dr sin r cos r u cos r cos r uj dr sin r ' 
. d sinr. 1 , d cosr, „ ik —ik 

-<»h ~ ) - h- + — /2 + /2 = , 

dr cosr w ar smr cosr uj cosr 

that is equivalent to the identity = and the equation for f 2 : 

k 2 



d 2 „ . cos r sin r . d „ , 9 1 

/ 2 + ^- h + (w - 1 

ar z smr cosr ar 



sin 2 r cos 2 r 



)/2 = 



(6.5) 
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The latter can be simplified: 

Mr) = — E(r) , |f - ^ + ( w ' = ; (6.6) 

sin r ar z sin r cos r ar cos^ r 

two concomitant functions are given by 

h(r) = — l — E(r) , / 3 = - — ^E(r) . 

u cosrsmr cu smr ar 



Turning back to eq. (|6.6p . first let us consider a particular case when k 2 

sin r d cos r d _ , 9 . 1 . _ 

cos r ar sm r ar cos z r 

from whence it follows 

cosr^ (^L± )E = k 2 E, 
sinr Or sinr dr 



so that 



cosr d d dr cosr 



sin r dr dx ' dx sin r ' 

-dx = d log cosr , x = log(C cos -1 r) , C = const 



and eq. (|6.7p reads 



d 2 „ 
—2 E = k 2 E , 

ax 



which has two solutions 

k 2 = u? , E± = e^ kx = E (cos r) ±k . (6.7) 



Among solutions 



(cosr)+ fe , re[0,£] (6.8) 

(cosr) /t 2 



at > the second one must be rejected because it turns to infinity as r — > 7r/2; when for < 
the first one must be rejected by analogous reason. Thus, physical solutions are 

k 2 = lo 2 , k > , E = E cos k r e -*M- fc2 ) ; 

k 2 = uj 2 , k < , E = E cos~ k r e^"*"**) . (6.9) 

Because at r / 0,tt/2, the values z = —tt and z = +tt determine one the same point in 
spherical space 53, solutions (|6.9p represent continuous function in S3 only if k takes on integer 
values: 

k = ±n, n = 1,2,3,... ; (6.10) 



9 



up c 
k = ± — = n , uj = — n , n = l,2,3,... (6.11) 
c p 



or in usual units 



Turning again to eq. (|6.6p let us construct one more special solution. To this end, considering 
approximate solution in the vicinity of the point r = 0: E = sin" 4 r, we get ^4 = 0, +2. In the 
same manner nearby the point r = ir/2 an approximate solution is E = cos B r , B 2 = k 2 . Let us 
demonstrate that there exist values k such that an exact solution can be constructed as follows 
(the case B = was considered above): 

E = sin 2 r cos B r . (6-12) 

With this substitution, eq. (|6.6p gives 

2 cos 4 r — 2(B + 1) sin 2 r cos 2 r — 3B sin 2 r cos 2 r + B(B — 1) sin 4 r 

2 2 "^2 "^2 2 

—2 cos r + B sin r — k sin r + sin r cos r = . 
With notation cos 2 r = x it can be written as 

2x 2 + (x- x 2 ) [-5B - 2 + uj 2 } + (£ 2 - B)(l - 2x + x 2 ) - 2x + 5(1 - a?) - /c 2 (l - x) = , 
or 

x 2 (4 + AB - uj 2 + -B 2 ) + x (-4 - AB + w 2 - B 2 ) + x° (£ 2 - k 2 ) = . 
The latter is satisfied if 

B 2 = k 2 , (B + 2) 2 -uj 2 = , 

that is 

B = -2 + uj, -2-uj , 
= ±5 , E = sin 2 r cos 5 r . (6.13) 

Thus, the corresponding solutions of this type are 

E = sin 2 r cos B r e -*(«"*-**) . (6.14) 

Solutions with negative must be rejected because they give infinite electromagnetic field at 
the point r = ir/2. Besides, periodicity requirement on z leads to k = ±1, ±2, ±3, .... 
Therefore, the wave propagating in the positive direction is given by 

B = +k = +l,+2,+3,..., 
k = -2 + u, uj = 2 + k = 3,4,5, ... , 

E = sin 2 r cos k r e -*M-M . (6.15) 

In turn, the wave propagating in the negative direction is given by 

B = -k = +l,+2,+3,..., 
-k = -2 + uj , uj = 2 - k = +3, +4, +5, ... , 

E = sin 2 r cos~ k r e -*M~M . (6.16) 
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Turning to the general equation (|6.6p . one may try to construct all other solutions of that 
type m = on the base of the following substitution: 



E(r) = sin A r cos B r F(r 



eq. (I6.6h gives (let cos 2 r = x) 



d 2 d 
Axil - x)—rF + 4[1 - 3x + B(l -x)}— + 
dx A dx 



+ 



„s 2x . , 1 — x 

- 2 + 5B) + - + B(B - 1 

1 — x x 



B k 2 2 
+ + uj 1 



1 — X X X 



F = 



Requiring k 2 = B 2 , for F we get the equation 



4x(l - x)^F + 4[(1 + B) - (3 + B)x] ^ - [ (B + 2) 2 - uj 2 ] F = 
which is of hyper geometric type 

z{l - z) F + [7 - (a + (3 + l)z] F' - a/3 F = , 
B + 2-LU „ S + 2 + w 



k = ±B , 7 = 1 + 5, a 



(6.17) 



(6.18) 



2 ' ' 2 
Thus, the general solution of the type m = takes the form 

E = sin 2 r cos B r F{a, /3, 7, cos 2 r) e ~ l M-M . 

We are to separate single- valued and continuous functions in S3. 
For a wave propagating in the positive direction z: 

k>0, k = +1, +2, +3, ... ; 

the function E(r) is finite at r = n/2 only if B = +k, besides polynomial solutions arise only if 

k + 2 - co 



(6.19) 



a 



-n =0,-1,-2,... =$> 
uj = k + 2(n + 1) = N . 



(6.20) 



For a wave propagating in the positive direction z: 

k<0, -k = 1,2,3,... ; 



(6.21) 



the function E(r) is finite at r = n/2 only if i? = — k, additionally one must obtain polynomials 
which leads to 



-k + 2 - uj 



a 



-n =0,-1,-2,... 



UJ 



-k + 2(n + 1) = N . 



(6.22) 



All constructed solutions of the Maxwell equations are finite, single- valued, and continuous 
functions in spherical Riemann space S3. 
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7 Maxwell solutions at k = 



Turning to eqs. (|6.3p at fc = 0: 

, d cos r sin r , „ im „ 

h- + /l + — /2 = , 

ar smr cosr smr 

„ im 

h — : — J"3 j 

a; smr 

I d sin r 

J 2 — H J 73 , 

cj ar cos r 

, d cosr, im „ . 

-w/3 + h- + — /2 - — /i = . (7.1) 
ar sin r sm r 

With the use of second and third from first and fourth we get an identity = and the 
following equation for f%: 

d 2 „ .cosr sinr. d „ , 1 k 2 . „ , 

7-2/3+ — / 3 +w 2 -l 5 — /3 = 0. 7.2 

dr z sm r cos r ar cos z r sm r 



which gives 



1 , d 2 E 1 dE 



771 



2 



h{r) = E(r) , — + — — h (to - — 2 )E = . (7.3) 

cos r ar z sm r cos r ar sin r 

First, consider a particular case m 2 = uj 2 : 

d 2 E 1 dE „ 1 , . . 

-T2 + + m 1 - s = ; 7 - 4 

dr z sm r cos r dr sin r 

with the solutions 

Physical solutions are 

m 2 = uj 2 , m > , E = E sin m r e^^^ ; 

m 2 = uj 2 , m < , E = E sin~ m r e -*(«"*-"»*) . (7.6) 

Performing analysis like in previous Section, we easily construct solutions: 

B = +m = +l,+2,+3, 
m = — 2 + uj , u> = 2 + m = 3,4, 5, ... , 

F 02 = cos 2 r sin m r e -*(<"*-™*) . (7.7) 



and 



B = -m = +l,+2,+3,..., 
-m = —2 + uj , ui = 2 — m = +3, +4, +5, ... , 

F 02 = cos 2 r sin"™ r e -*(<"*-"") . (7.8) 
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All possible solutions of eq. (|7.3p can be constructed on the base of a substitution: 

E(r) = cos 2 r sin 5 r F(r) . (7-9) 

and further (let sin 2 r = x) we get 

d 2 d 
m 2 = B 2 , 4x(l - x)-^F + 4[(1 + B) - (3 + B)x] — - [ (B + 2) 2 - lo 2 } F = , 

what is of hypergeometric type 

z(l - z) F + [7 - (a + p + 1>] F' — a/3 F = , 
k = ±B, 7 = 1 + 5, a = - g + 2 ~ u; ; p 



2 7 ' 2 
Thus, the Maxwell equations solutions of the type k = is given by 

E = cos 2 r sin 5 r F(a, /3, 7, sin 2 r) e^M-"^) . 

We are to separate physical waves. 
m > , B = +m , 

m + 2 — w 



m < 



a - 
B 



a 



= —m, 

-m + 2 — uj 



-n 



-n 



to = m + 2(n + 1) = N 



(7.10) 



(7.11) 



-m + 2(n + 1) = N 



(7.12) 



(7.13) 



8 Radial system at arbitrary m, /c, general solutions 



Now let us solve radial equations in general case (]6.3p . The first equation reduces to the identity 
= when taking into account three remaining: 



ik 



h 



im 



h 



cos r sin r 
d sin r ik 

-uh = {-, ) h h 

ar cos r cos r 

d cosr im 
-w/ 3 = -(3- + ttt- ) /2 + ITTt/i ! 



and the first equation takes the form 
d cosr sinr ik 

\T + ~ )( — J 2 

ar sm r cos r cos r 



' dr sin r 



h) + 



sm r 



smr 



smr 



d sinr ik 

{-j ) J3 h 

ar cos r cos r 



+- 



ik 



cosr 



d cosr 

-("J" + — ) 72 + 7i 

ar sm r sm r 



+ 
. 



what is the identity = 0. The system (J8JJ) is simplified: 



h 



1 



smr 



1 



cos r 



(8.1) 
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so that 

k F 2 — m F3 



u fx = i 



sin r cos r 



F 2 1 dF 3 ik 

-oj - — = fi , 

smr cosr dr cosr 

-.^ = -J-fi + 4^/i. (8.2) 
cosr smr dr smr 



Excluding f\, we arrive at 



. lu d km . _ 1 . 9 k 2 . ^ 

J- + 2- ^3 + — ^ J- ^2 = 0, 

cos r ar sin r cos z r sin r cos z r 

. a; d km N „ 1,9 m 2 , „ , 

(— j — F 2 + (-w 2 + ~ ; 2 ) ^3 = 0. 8.3 

sm r ar cos r sin r cos r sin r 

With the use of a new variable y = (1 — cos 2r) /2 the system reads 

(2 ^ ^ m ) i 7 + ( ^ + m2 ) i 7 

dy y{l-y) l-y y(l - y) 

d km , „ . u 2 A; 2 , „ . 

< 2 ^ + ^»^ + < + 7-ra )F2 = ' <8 - 4) 

Instead of i*2 , -F3 let us introduce new functions by means of linear transformation with unit 
determinant aN — 0M = 1: 

F 2 = a(y) G 2 + (3(y) G 3 , 
F 3 = M(y) G 2 + N(y) G 3 , (8.5) 

and inverse given by 

G 2 = N(y) F 2 - p(y) F 3 , 
G 3 = -M(y) F 2 + a(y) F 3 . (8.6) 

Combining eqs. (|8.4p we get 

d km . ^ dN ^ diV „ w 2 m 2 , „ 

AT (2w- r F 2 + 2w— F 2 - 2uo—F 2 + iV -- + — F 3 - 

dy y{± -y) dy dy l-y y{\ - y) 

dy y(l - y) dy dy y y(l - y) 

, r ,d km . „ dM „ dM „ „„ w 2 m 2 , „ 

-M (2w- r F 2 - 2w— F 2 + 2w— F 2 - M(— + — F 3 + 

dy 2/(1 - y) dy dy l-y y(l - y) 

d /cm , da _ da _ . a; 2 k 2 . „ , . 

+a 2w— + — F 3 + 2w— F 3 - 2uo—F 3 + a + A F 2 = , 8.7 

dy y(l - y) dy dy y y(l - y) 
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from whence it follows that 



d „ km ^ diV _ . w 2 m 2 . _ 

2a; G 2 - iV — F 2 - 2u—F 2 + N -- + — F 3 - 

dy y{± -y) dy 1 - y y(l - y) 

km „ dd „ „ . A; 2 

2/(1 -y) ay y y(l - y) 



„ d „ „ r fem _ „ dM _, , , . to m 

2lo— G 3 + M — F 2 + 2lo—F 2 - M -- + — 

dy y{± -y) dy 1 - y y(l - y) 



F 3 + 



km „ da „ . a; fe . 
+a — F 3 - 2lo—F 3 + a + F 2 = 

-y) ay y y(i-y) 



Instead of F 2 , F 3 we substitute their expression through G 2 , G 3 according to 

km 



2lo— h 

dy 



- (Na + 0M) -— — -2^a + NM ~^ + r f 

y(i-y) «y y(i-y) 



+ 



d/3 



oj 



\i-y)-k 2 



dy y(l - y) 



G 2 + 



+ 



-2NP 



km 



dN 



y(i-y) y(i-y) 



d/3 Ar ,2 u; 2 (l- y )-k 2 



+2uj^-N - p 
dy 



y(i - y) 



G 3 = 0, 



9 

2uj— h 

dy 



» T x km dM „ — u; 2 y + m 2 

M/3 + aAT — + 2a;— P-NM — — ■ 

2/(1 - y) dy y(l - y) 



da 



'(l-y)-k 2 



-2oj—N + Pa 

dy y(l - y) 



G 3 + 



+ 



km dM Q —u 2 y + to 2 
2Ma — r + 2w— a - M — — - 

y(i-y) ay y(i-y) 



da , , 
-2u)—M + a' 
dy 



u 2 (l-y) - k 2 



2/(1 - y) 

Let us assume that the transformation used is an orthogonal one: 

a G 2 + P G 3 = cos A G 2 + sin A G 3 , 
M G 2 + N G 3 = - sin A G 2 + cosAG 3 , 



G 2 = 0, 
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then 



and 



- 2uj— a + 2uj-^-M = -2u[(cos^)'cos^ + (sin A)' sin A] = , 
ay dy 

dN dd 

-2uj— (3 + 2u-t-N = 2uj[ - (cos AY sin A + (sin A)' cos A = + 2u A' , 
dy dy 

2oj— f3 - 2uj—N = 2w[-(sin J 4) / sinA - (cos A)' cos A] = , 
dy dy 

2oj^- a - 2u-^-M = 2w[-(sin,4) / cos,4 + (cos A)' sin A] = - 2u A' . 
dy dy 



Na + j3M = cos 2A , 2N(3 = sin 2A , 2Ma = - sin 2A , 
a(5 = sin A cos A = - sin 2A , NM = - sin A cos A = — - sin 2^4 , 
jV 2 = cos 2 A , /3 2 = sin 2 ,4 , M 2 = sin 2 ,4 , a 2 = cos 2 ,4 , 



Therefore, eqs. (|8.9p and (|8.10p take the form 
dG 2 



2w- 



c/y 



nA km 1 . -w 2 y + m 2 + w 2 (l-y)-/c 2 
cos2A^^ 7 H — siniA 



2/(1-2/) 2 



2/(1 - 2/) 



G 2 + 



+2^71' - sin 2A 



km 



2/(1 - 2/) 



+ cos z A 



2 , -w 2 2/ + ^ 2 



2/(1 - 2/) 



sin 2 74 



w 2 (l-y)-/c 2 
2/(1 - 2/) 



G 3 =0, 



2uj—± + 
dy 

-2uA' - sin 2 A 



cos 2,4 



km 



+ - sin 2^4 



-uj 2 y + m 2 + w 2 (l — y) — k 2 



2/(1-2/) 2 

/cm 



2/(1 - y) 



G 3 + 



2/(1 - 2/) 



.-, , —oj z y + m z 9 , oj 2 (1 — y) — k 2 

sin 2 A ^ — + cos 2 yl v 



2/(1 - y) 



2/(1 -y) 



G 2 = 0, 



Supposing that the used linear transformation does not depend on coordinate y, we get more 
simple expressions: 



2lo- 



dG 2 
dy 



nA km 1 . -uj 2 y + m 2 +uj 2 (l-y)-k 2 
cos 2A^-_ r H — sm 2A 



2/(1-2/) 2 



2/(1 - 2/) 



G 2 + 



-2/cmsin2^4 + (1 + cos2A)[-ui 2 y + m 2 ] - (1 - cos 2,4) [oj 2 (1 - y) - k 2 ] 

+ M^V) G3 = ' (8 ' 12) 



dG 3 
dy 



cos 2A 



/cm 1 — a; 2 ?/ + m? + u; 2 (l — y) — k 2 1 

— H — sin 2 A ; 

2/(1-2/) 2 1/(1-1/ 



G 3 + 



—2km sin - (1 - cos2^)[-w 2 y + m 2 ] + (1 + cos2^)[w 2 (l - y) - k 2 ] 

+ m^y) ° 2 = 



us) 
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Let 



cos 2 A = , 



2A = - 
2 



sin 2 A = 1 



then eqs. (l87^+|8+3l) read 



d -u} 2 y + uj 2 (1 - y) + m 2 - k 2 \ -J 2 + (m - k) 2 

2U Ty 2W^V) ) ° 2+ 2y(l-y) ° S = ° 



^ <i — w 2 y + w 2 (l — y) + m 2 — k 2 
dy 2y(l - y) 



u; 2 -(m + k) 2 
° 3+ 2y(l-y) G2 "°- 



(8.14) 



It is remarkable that in the system produced the singularities are located at the points 
y = 0, l,o o) on ly. 

From (|8.14p and excluding G3 one straightforwardly gets the equation for G 2 : 



G, 



-2a;- 



2y{\-y) dG 2 -uj 2 y + uj 2 {1 - y) + m 2 - k 2 



-to 2 + (m — k) 2 dy 



+ 



%(i-y) 



d 2 Go . , s dG 2 / 

+ 4(1-22/)— - + [2lo + u; 



a; 2 + (m — k) 2 

- m 2 — /c 2 
+ 



G 2 



G 2 = 



dy 2 " 3/ ' V"" ~ viX ~ y) ' 1-2/ 

With the substitutionG 2 = y A (l - y) B G(y), eq. (|8.15p takes the form 

4y(l - y) G" + 4 [A(l - y) - By + A(l - y) - By + (1 - 2y) ] G ' + 



+ 



4A(v4 - 1)- + 45(5 - 1)^— 

y 1 - y 



44(A - 1) - 45(5 - 1) - 8AB+ 



+4(-2A - 2B + - + 

y 1 - y 



.1 



+ 2w + or - m (- + 



y 1-2/ 



) + 



A' 2 



1-2/ 



G = 



Requiring 



4A(,4 - 1) + 4A - m 2 = 
45(5 - 1) + 45 - k 2 = 



1 , 

A = ± — m 

2 1 

B = ± - I fe I : 
2 1 1 



we arrive at 



y(l - y) G" + [2A + 1 - 2{A + B + l)y ] G ' - 
A + B)(A + 5 + l)-|(| + l)] G = 
what is of hypergeometric type 

7 = 2A + 1 , a + /3 = 2,4 + 25 + 1 , 
a/?=(A + 5)(A + 5 + l)-|(| + l), 



(8.15) 



(8.16) 



(8.17) 
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that is 



a = A + B-~, @ = A + B + 1 + ^ , 7 = 2A + 1 
The functions are finite on the sphere S3 only if 



(8.18) 



1 , 

A = +— m 

2 1 



B = + \\k 



a 



-n = 0,-1,-2,... ; 



which leads to the frequency spectrum in the form : 

co = 2(n + A + B) = 2n+ \ m \ + \ k | ; 



(8.19) 



the parameters m and k are allowed to be integer only : m , k € {0, ±1, ±2, ... }. The function 
G 2 {y) is 



G 2 (y) = M 2 y H/2 (1 - y) |fc|/2 F(-n, n+l+\m\ + \k\, \m\ +1; y) . 
In eqs. (|8.14|) we might exclude GV G3: 

^ , s 2y(l - y) dG 3 -LU 2 y + u 2 (l-y)+m 2 -k 2 ^ 

^2{y) = —lU) — —TT^~ 5 7 ! TTo ^3, 



cj 2 — (m + A;) 2 dy 
4y(l-y)^ + 4(l-2y)^+(-2^ 2 



u; 2 — (m + k) 2 



m 



+ 



m 



k 2 



2/(1-2/) i-y 

With the use of substitution G3 = y^(l — y) B F(y), the equation for G3 reduces to 
4y(l - y) F" + 4 [A(l - y) - By + A(l - y) - By + (1 - 2y) ] F ' + 
+ 



4AM - 1)- + 4£(£ - 1)— 4A(A - 1) - AB(B - 1) - 8AB+ 

2/ 1-2/ 



+4(-2vl-2 J B + - + -?-)-2^ + ^-m 2 (i + -^) + ^-^ 

y 1 -y y l-y 1-2/ 

Requiring 

1 

AA(A - 1) + 4A - m 2 = = 
4£(B - 1) + 45 - k 2 = = 

we arrive at a hypergeometric type equation 

y(l - y) F" + [2A + 1 - 2{A + B + l)y ] F ' - 

F = 



F = 



A = H — I m 
2 1 

B = + - I Jfe I : 
2 1 



(A + J B)(A + B + 1)-|(|-1) 



a = A + B + 1 



c = 2A + 1 



(8.20) 



G 3 = (8.21) 



(8.22) 



(8.23) 
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Further we get 



a = A + B + l-~ = -N, AT = 0,1,2,..., 



lo = 2(A + B + 1 + N) =| m \ + | k \ +2(1 + N) , N + l = n , 

G 3 = M 3 yH/2 (1 _ y) |fc|/2 F( _ n + i, „+ | m | + [ |, | m | +1; y) ; (8.24) 

compare with (I8,20p , 

It remains to find a relative factor in two functions G 2 and G3: 

G 2 = M 2 y H/2 (1 - y) |fc|/2 F(-n, n + 1+ | m \ + \ k |, | m | +1; y) , 

G 3 = M 3 y H/2 (1 - y) |fc|/2 F{-n + 1, n+ | m | + | k |, | m | +1; y) , (8.25) 

and the relationship (see (|8.15p ) 

G 3 [(m - fc) 2 - oo 2 ] = -Aw y(l - y) ^ + [m 2 - k 2 + uj 2 (1 - 2y)\ G 2 . 

ay 

must hold. Using the expressions for G 2 G 3 we get 

M 3 

(m-k — 10) (m-k + 10) — F 3 (y) = 

M 2 

= -Auj [ J-^i (1 - y) F 2 (y) - IA1 y F 2 (y) + 
+y(l " y) ^ ^2(y) ] + [m 2 - /c 2 + w 2 (l - 2y)] F 2 (y) . (8.26) 

It is sufficient to consider this equation in the point y = only that results in 

M 3 

— (to + m — k) (to — m + k) = (u— \ m \ —k) (uj— \ m \ +k) , 

M 2 

and therefore 

M 2 = M {uj + m- k){u -m + k) , 
M 3 = -M (lo- I m I -k)(uj- \ m\ +k) . (8.27) 

Depending on the sign of m it may be rewritten in a simpler form: 

m > , M 2 = M(uj -k + m) , M 3 = -M(co - k - m) ; 

m < , M 2 = M(co + k-m) , M 3 = -M(w + fc + m) ; 

m = , M 2 = M , M 3 = -M ; (8.28) 

M stands for a numerical constant. 
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9 Maxwell solutions in elliptical model 

Let us consider the problem of Maxwell solutions in elliptical space S 3 . This space S 3 is 
a space of constant positive curvature also and differs from the spherical model in topological 
properties only: S3 is 1-connected, S 3 ia a 2-connected. The question is on the role of these 
differences for electromagnetic field solutions. 

To obtain explicit realizations for two models it is convenient to use relations known in the 
theory of unitary and orthogonal groups. To each point in S3 there exists corresponding element 
in unitary group SU(2): 

B = cr° n — i a k n k , det B = +1 . 

In turn, to each point in elliptic space S 3 there exists corresponding element in SO (3) 
parameterized by Gibbs 3-vector [23] 

0(c) =1 + 2 -A— > , (c x ) kl = - e klj Cj . (9.1) 

note that two infinite length vectors represent one the same point in S 3 : 

c+ = +00 c , c - = -00 co, c 2 = 1 , 0(c ± 00 ) = / + 2 (c x ) 2 . 
Mapping 2 — > 1 from SU(2) to SO(3) is given by 

{+n a ; -n a \ -> c = — . 

n 

Cylindric coordinates (p, <fi, z) in elliptic space can be defined by the relations 

tan p tan p . 

c\ = cos z , C2 = sin z , C3 = tan z , 

cos p cos 2 

G , p G [0, tt/2] , G [-7T, +7r] , Z G [-7T/2, +7T/2] . (9.2) 

Additionally, we must define such an identification rule on the boundary of the region G, which 
agrees with the identification rule for vectors cj, and c^. To this end, it is convenient to divide 
the region G into three parts: 

Gi = G(p + 0,vr/2) , G 2 = G(p = 0) , G 3 = G(p = tt/2) . 

For the region G\ identification is given by (for more detail see [...]) 

Fig. 1 Gi 



z 





1' B 2 


+7T/2.4 


1" 


c 






G' 










D 


















2' A' 1 


-it/2B' 


2" 
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each pair (A, A') , (B,B') and so on represents one the same point in elliptical model S' 3 ; also 
(1, 1', 1") and (2, 2', 2") correspond to one respective point in S' 3 . 

Now we should find what of above constructed Maxwell solutions in case of spherical model 
S3 will be single-valued ones when considering elliptical model S' 3 (here we examine only points 
parameterized by the region Gi). Evidently, it is sufficient to examine the behavior of the vector 
/ = e tm ^ e lkz . From the relations 

/(C) = /(C), f(D) = f(D'), ... 

no additional restrictions arise besides that M and K to be integer. Equations f(l) = f(l') = 
/(I") give 

g — ik(ir/2) _ g—imn e +ifc(-7r/2) _ e +imir e +ik(n/2) 

whence it follows that 

i2rmz -1 i(k—m)ir i i(k+m)n -i . 

O J- ^ O J- 2 C- -L « 

and therefore, (k — m) and (k + m) must be even. The same results follows from consideration 
of points (2, 2', 2"). 

Therefore, Maxwell solutions constructed above in spherical space will be single-valued so- 
lutions in elliptical space (in region G\ ) only if m and k are both even, or both odd; corre- 
spondingly, to parameter N in the expression for the frequency spectrum (see. (|8.19p ) given 
by 

uj = 2n+ I m I + I k \= N (9.3) 

takes on even values: N = 0, 2, 4, 6, ... 

Now let us consider the behavior of the above mentioned Maxwell solutions in the remaining 
regions C and G3. 

First, let us specify the case of Gi and consider the vicinity of the point P (see Fig. 2): 
p = (vr/2 tan a + z tan a) {z = —vr/2 + 5, p = 5 tan a} ; 



Fig. 2 


G2 (<? 


) is arbitrary) 






L P . 

+VT/2 




P 






P' 


-it/2 


+vr/2 
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and 



tan (5 tan a) cos (ft tan (8 tan a) sin < 

c l = ■ /r , \ : — F > c 2 — 



COS 5 



sin (5 tan a) sin (5 



C3 = : 



sin S 



sin(<5 tan a) sin 5 

so in the limit S — > (a / 0) we get 

G 2 , P, c= + oo(0, 0, -1) . 

This means that coordinate (ft is "mute" at the point (p = 0, z = — 7r/2, 0). 
For another point P' = (p = 0, z = +7r/2, 0) we get similar result: 



Go 



P' 



oo (0 , , -1) 



(9.4) 



(9.5) 



Thus, P and P' represent one the same point in elliptical space. Compare & wm k at these two 
points P and P'\ 



®umk ~ e im0 e lkz (sin r)H (cos r)l fe lF(A, 5, C; sin 2 r) ; 



1*1 . 



and 



0, if 0; 

e +ikn/2 F (A,B,C;0), if m = 0; 



0, if m ^ 



^mfc(P') ~ 



e -ifc7r/2 ^(^^(7.0), if m = ; 



(9.6) 



because k is even (at m = 0) we have the identity &u)mk{P) = *&umk{P')- ^ n remaining part of 
the region G2, points of elliptical model are parameterized according to 



(O;0;z^O,7r^2) 
c = (0 , , tan z) , G2 {(ft — "mute" variable) . 

and the function <& wm fc 



Go 



Yojmk 



0, if m + 0; 

e ikz f(A,B,C;0) , if m = , 



(9.7) 



(9.8) 



is single-valued and continuous in that part of G2 of the elliptic space. 
Finally, let us consider the region G 3 : 



Fig. 3 



G3 {(ft is arbitrary) 



B 




■ P 

+7T/2 


B' 




\ 


/ a 




-it/2 


+7T/2 
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In the vicinity of B we have 

7T 7T 7T 7T 

p = [ -z tan a + (- - - tan a) } -> { z = (-- + 6) , p = (- - 5 tan a) } ; 
that results in 

cos(<5 tan a) cos(/> cos(£ tan a) sine/) cos<5 

1 sin((5 tan a) sin 6 ' 2 sin(<5 tan a) sin 5 ' 3 sin 5 ' 

from whence in the limit 5 — > (as a 7^ 7r/2 ) it follows 

-B : c = 00 (cos0, sine/;, 0) = 00 (cos</>, sin^, 0) . 

tan a 

Analogously, in the vicinity of B' we get 

B' : c = 00 (cosc/>, sin</>, 0) . 
Comparing function at points B and B': 



, if k ^ 0; 

e +im7r/2 i?( y 4 5 5 5 C;l) , if k = 



<5>. mk {B')~{ l-JJt(A,B,C-,l), if = ; (9 ' 9) 



because m is even, the equality ^uimk(B) = $>umik(B') holds. In remaining part of G, 



3- 



G 3 : (vr/2, z + -ir/2, +tt/2) => 
00 

(cos0, sin^>, 0) = 00 (cos(/>, sin</>, 0) . (9.10) 



cos z 

z is "mute" coordinate, therefore the solutions 

^f0,iffc^0, . , 

3 \ e imz F(4,5,C;1) , if fc = , 1 ; 

represent single-valued and continuous functions in that part of elliptic space. 

Thus, all functions $uimk(p, (j>, z) at uj = N = 0, 2, 4, ... are single-valued and continuous in 
elliptic space, and they represent physical solutions for Maxwell equation in this space, whereas 
all remaining functions <&ojrnk{p-,4>i z ) > N = 1, 3, ... should be rejected as non-single-valued 
and discontinuous in S 3 space. 



10 Cylindric coordinate and tetrad in Lobachevsky space #3, 



Let us consider Maxwell equations in cylindric coordinate [22J of hyperbolic Lobachevsky space 
n\ = sinh r cos <j) , ^2 = sinh r sin <j) , 713 = cosh r sinh 2 , 724 = cosh r cosh z ; 



i)(y) 



(iS- 2 = dt 2 - dr 2 - sinh 2 r 

10 

10 

sinner 

cosher 



e(a)/9(tf) 



cosh 2 r (fe 2 , x a = (i, r, 0, z) , 

10 

0-10 

— sinh r 

— cosh r 



(10.1) 
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where (r, (ft, z) run within 

re[0,+oo), 0G[O,27r], z€(-oo,+oo) 
Christoffel symbols are 



"/3<r — , 



oo 



0, 



0j 



0, 



jk 





— sinh r cosh r 
— sinh r cosh r 



1 jfc 



'- 1 jk 







sinh r 
cosh r 



cosh r 



Derivatives of tetrad vectors read 



e (0)/3;a - -g^- ~ r a/3 e (0)(T - , 



sinh r 



cosh r 







cosh r ri 
sinh r 





sinh r 





5e (l)/3 



e(l)/J;a = -q^T - r a/3 e (l)r 
















sinh r cosh r 






sinh r cosh r 



de 



(2)/3 r 



TV 

e (2)/3;a - - 1 a/3 e (2)<^ 



<9e 



2(3)/J;a 



(3)0 



-r^ e 



a/3 e (3> 


















cosh r 



































sinh r 





















Ricci rotation coefficients are 

7011 = 7021 = 7031 = , 7012 = 7022 = 7032 = , 7oi3 = 7023 = 7033 = 



and 



7231 = , 7311 = , 7121 = , 

cosh r 

7232 = , 7312 = , 7122 = -T~, , 

smnr 

sinh r 

7233 = , 7313 = — , 7123 = . 

coshr 
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With relations 



d, 



1 



7 (0) 



dt, 



'(2)^P °(' 2 ) ginhr 



i 



9, 



(3)^7* ( 3 ) coshr" 



Vo = (7010,7020, 7030 ) = , 
V2 = (7010 , 7022 , 7032) = , 

Po = (7230,7310,7120) = , 



d z , 

1 r 

vi = (7011,7021,7031) = , 

V3 = (7013,7023,7033) = , 

Pi = (7231,7311,7121) = , 
cosh r , 



P2 = (7232,7312,7122) = (0,0, 
P3 = (7233 , 7313 , 7123) = (0 



sinhr 
sinhr 



) 



Maxwell equation (|16.3p reads 



-id+ + a 1 d r + a 2 



1 



sinhr 



1 



a. ' 



cosh r 



d z + a 2 S : 



coshr 
' sinh r 



a 3 S: 



coshr 



sinhr 

> 

' cosh r 



0) 





E + icB 



0(10.3) 



11 Separation of variables in space H3 



Maxwell matrix operator in f)10.3|) commutes with three operators: idt, id^, id z ; therefore 
solutions can be constructed on the base of the following substitution 





E + icB 





fi(r) 
/ 2 (r) 
/a(r) 



and eq. f)10.3j) takes the form 



-oj 



1 d 
a 1 — + 



im 



dr sinh r 



ar + 



ik 



coshr 



a 3 + 



coshr 2 sinhr o 

— — a 63 — a°b 2 

smn r cosh r 



After simple calculations we get the radial system in the form 

d ^ coshr ^ sinhr ^ im ^ ik 
dr sinhr coshr sinhr coshr 



ik 



im 



h 

h 



cosh r sinh r 
d sinhr ifc 
- (3- + — — ) h + ttzt—Zi 



'<ir ' coshr' 1,0 coshr 
„ . d coshr im 

-w/ 3 + + -r-r-) /2 - — T— 

dr smn r smn r 



/i 



0. 

o, 
o, 

= . 





fi(r) 
hi?) 
hi?) 



11.11 



0. (11.2) 



111-3) 
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12 Solutions at m = 0, in space H% 

First, let m = - then 

, d coshr sinhr. „ ik 
ar smnr cosnr cosnr 

f - ~ ifc f 

Jl — 7— J2 , 

uj cosnr 
, d sinhr. „ ik 

-vf2 ~ h- + — r- /s + — r-/i = , 
ar cosn r cosh r 

1 , d coshr N „ 

h = ~ h- + -!- / 2 • 12.1 
a; ar sinn r 

With the help of 2nd and 4th equation from the 1st and 3rd it follows the identity = and 
equation for fi'- 

„ , , 1 , . d 2 £ 1 dE , 9 k 2 . . 

f 2 (r) = —- E(r) , -fY-^r, T~ + ^ Z2 ) = ' ( 12 - 2 ) 

sinn r dr z sinn r cosn r ar cosh r 

besides 

/i(r) = — h \ r £?(r) , / 3 = - £i*(r) • 

a; cosnr smnr u smnr ar 

Eq. (|12.2p can be resolved very easily when k 2 = uj 2 : 

sinh r d cosh r d „ 9 . 1 . „ 

^^^-1—^ E + fc 1 ^ = ; 

cosn r ar sinn r dr cosh r 

from whence it follows 

c^d c^d 
sinh r dr sinh r dr 

With the help of a new variable 

cosh r d d dr cosh r 



sinh r dr dx ' dx sinh r ' 

dx = d log cosh r , x = log(C cosh r) , C = const . 

we arrive at 

d 2 

— ? £ = —k 2 E , E = e lkx = const (cosh r) lk , k = ± to ; 
dx z 

Two constructed solutions are conjugated: 

(cosh r) +ik = ( e lo s cosh r )+ ik = cos ( k log cosh r)+i sin( A; log cosh r ) , 
(cosh r)~ ifc = [ e lo s cosh r ]~ ik = cos ( k log cosh r ] - i sin( log cosh r ) , 

therefore one can separate two independent real ones: 

E + (r) = cos [ ko log cosh r ] , E^(r) = sin [ ko log cosh r ] . (12-4) 
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In the limit of vanishing curvature they reduces to the known ones : 

r — ► , E + (r) = cos ( ko log cosh r ) — > +1 ; 
r — »■ , E-(r) = sin ( &q log cosh r ) — > const r 2 ; (12.5) 

these satisfies to the equation in the flat space-time 

4--4~ e (t)=0, =► £(r)~l, r 2 . 
ar r dr 

In contrast to the flat space, in Lobachevsky model the waves (I12.4|) are both oscillating at 
infinity. 

Let us show that (112.21) has another simple exact solution in the form 

E = sinh 2 r cosh 5 r . (12.6) 

Indeed, substitution of (112.60 into (|12,2p leads to 

2 cosh 4 r + 2(B + 1) sinh 2 r cosh 2 r + 3B sinh 2 r cosh 2 r + B(B — 1) sinh 4 r — 
—2 cosh 2 r — B sinh 2 r — k 2 sinh 2 r + u 2 sinh 2 r cosh 2 r = . 

With notation cosh 2 r = x it reads 

x 2 (4 + 45 + u 2 + 5 2 ) + x (-4 - 45 - u 2 - B 2 ) + x° (B 2 + k 2 ) = . 

The latter is satisfied if 

B 2 = -k 2 , (B + 2) 2 + to 2 = , (12.7) 

that is 

B = -2 + iw, -2-iw, fc = ±i 5 = { T !^ + ^! ' (12.8) 

\ T (2t - w) , 

Corresponding solutions look as 

r, z) = £ sinh 2 r cosh B r e" !(wt " fe) . (12.9) 

their real and imaginary parts are given by 

B = —2 + iu, k = iB = —2i — u, 

E(t, r.z) = E sinh 2 r cosh- 2+i w r e *(- 2 **- w *- wt ) = Eq tanh 2 r cosh^ r e 2* e i(-*«-w*) = 
= -Eo tanh 2 r e 2z [cos(u; log cosh r) + i sin(w log cosh r)\ [cos(— uz — ut) + i sin(— uz — ut)] = 
= Eq tanh 2 r e 2z cos(o; log cosh r — ujz — ut) + iEq tanh 2 r e 2z sin(u log cosh r — ujz — ujt) 

B = -2 + iuj, k = -iB = 2i+u, 

E(t, r.z) = E sinh 2 r cosh- 2+ * w r e *(2*i+«*-wt) = E Q tanh 2 r cosh^ r e -^ e ^~^) = 
= Eq tanh 2 r e~ 2z [cos(u; log cosh r) + i sin(w log cosh r)] [cos(uiz — ujt) + i sin(u;2! — ut)] = 
= Eq tanh 2 r e~ 2z cos(u log cosh r + uz — ut) + iEq tanh 2 r e~ 2z sin(u; log cosh r + uz — ut) 
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B = —2 — ioj, k = iB = —2i + oj, 

E(t, r.z) = E sinh 2 r cosh" 2 "* w r e <(-2*i+«*-wt) = Eq tanh 2 r cosh -^ r ^z^z-^t) = 
= Eq tanh 2 r e 2z [cos(— oj log cosh r) + i sin(— oj log cosh r)] [cos(a;z — u>i) + i sin(wz — ujt)\ = 
= Eq tanh 2 r e 2z cos(— oj log cosh r + ojz — ojt) + iEq tanh 2 r e 2z sin(— oj log cosh r + ojz — cot) 

B = -2 - iu, k = -iB = 2i - oj,] 

E(t,r.z) = Eq sinh 2 r cosh" 2 "* w r j(™-"*-<*t) = Eq tanh 2 r cosh"^ r e" 2 ^^^ 2 "^ = 
= Eq tanh 2 r e~ 2z [cos(— oj log cosh r) + i sin(— oj log cosh r)] [cos(— ojz — ojt) + i sin(— ojz — ojt)] = 
= Eq tanh 2 r e~ 2z cos(— oj log cosh r — ojz — ojt) + iEq tanh 2 r e~ 2z sin(— oj log cosh r — ojz — ojt) . 

Physical sense of these waves is not clear. More insight can be reached when constructing 
all possible solutions of eq. (|12.2p through the substitution 

E = sinh 2 r cosh 5 r F(r) . 

In that way, from (|12,2j) we arrive at an equation of hypergeometric type (if x = sinh r, k 2 + 
B 2 = 0) 

d 2 d 
Ax{\ - x)^F + 4[(1 + B) - (3 + B)x] — - [ (B + 2) 2 + to 2 ] F = , (12.10) 



with 



that is 



and 



1 + B, a + p = 2 + B, a p=^ + 2 ^ +UJ 



4 



B + 2-ioJ „ B + 2 + ioj 
B = ±ik, 7 = 1 + 5, a = , (3 = , (12.11) 



E(t, r, z) = sinh 2 r cosh B r F(a, 0, 7, cosh 2 r) e -*(<^-**) . (12.12) 

Evidently, above constructed solutions (|12.7p - (|12.9p can be obtained from the general 
relations (|12.1ip - (I12.12p . if one demands a = or (3 = 0. However, no physical ground exists 
to impose such (polynomial) restrictions in the case of Lobachevsky space. Instead, the complete 
electromagnetic basis should include waves spreading along z, with real parameters k. 



13 Solutions at k = in space H s 

Now let us turn to eqs. (jll.3p with k = 0: 



. d cosh r sinh r . „ im 

J" + —7- + ~T )/l + —7- h = , 
dr smnr cosnr smnr 

, im 

Jl — r-r— J3 , 

oj sinn r 

1 d sinh r 

h = {-7- H ; — J 73 , 

cj dr cosn r 

, <i coshr, im „ . 

-w/3 + h- + — r- h ~ -1-/1 = • 13.1 
dr sinn r smn r 
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The first and fourth equations gives 



, d coshr sinhr 
(-T + -^ + 



im im 
h + 



1 . d sinh r , „ . 



dr sinhr coshr u> sinhr sinhr u> dr coshr y 
. d cosh r , . 1 . d sinh r , „ , im im 
+ {-r + — rH (— (t- + — [-) /a) " irzrr: rriz^ h = . 



v <ir sinh r y v u> y dr ' cosh r 
they reduces respectively to the identity = and 



sinh r ui sinh r 



/s(r) 



1 



£(r) ; 



d 2 £ 



1 



dE 



+ — 



coshr v ' ' dr 2 sinhr coshr dr 
In variable y = — sinh 2 r it is rewritten as 



+ (u, 2 



ra 



sinh 2 r 



)E = . 



d 2 



d 



m 



%(1 - v)-^ £ - 4(1 - y)— £ + (u/ + — ) E = , 
dy 2 dy y 



that is solved in hypergeometric functions E = y a (l — y) b Y(y): 

d 2 

4y(l - y)- r5 Y + 4[1 + 2a - (2a + 26 + 1) 

dy 2 

-4(a + 6) 2 -a; 2 + (4a 2 - m 2 )- + 46(6- 1" ' 



dy 



Y + 



1-2/ 



y = o 



Requiring 

m = ±2a , 6 = 1 
we get an equation of hypergeometric type with 

7 = 1 + 2a, a + /3 = 2a + 26, 

that is 



a/3 



0, 



4(a + 6) 2 +w 2 



7 = 1 + 2a 



a = a + 6 =p 



/3 = a + 6± 



la; 



;i3.2) 



(13.3) 



(13.4) 



13.5) 



13.6) 



14 Solutions with arbitrary ra, k in space #3 



Now let us consider radial equations in general case fll 1 .3j) : the first equation in (jll.3p turns to 
be the identity = when three remaining hold: 



with substitutions 



ik 



h 



im 



/3 , 



cosh r sinh r 
d sinhr ik 

= (v- H 7— ) J3 7— Jl , 

dr cosh r cosn r 

, , d coshr im 

-uf 3 = -(— + — ^) / 2 + -Tzrzh ; 



v dr sinh r ' 



1 

sinhr 



sinh r ' 



1 



cosh r 
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one obtains 

k F 2 — m F 3 



- u fi = i 



sinhrcoshr ' 



-to F<2 = 1 dFs - ik U 
sinhr coshr dr coshr ' 

F3 _ 1 dF 2 | im 

coshr sinhr dr sinhr 



After excluding of f\ we get 



/ <jj d km . _ , u) 2 k 2 . „ 

cosn r dr smh r cosh r sum r smh r cosh r 

/ w d /cm , „ . uj 2 m 2 . „ , 

) F 2 ~ (— T- - — T-T-T2-) ^3 = 0. (14.2) 



sinhr dr coshrsinh 2 r coshr coshrsinh 2 

In variable cosh 2r — 1 = 2y, eqs. (|14.2|) take the form 

d km s ^ _|_ / ^ 2 _j_ m2 s ^ _ q 

dy y(i + y) i + y y(i + y) 

d km , , uj 2 k 2 , , 

( 2u 7T + 7TZ \ ) F * + + n— T ^ 2 = • 143 

^ + y) y y(i + y) 

Let us introduce new functions by means of a linear transformation (with unit determinant 
aN — (3M = 1 and numerical parameters) 

F 2 = a G 2 + G 3 , 
F 3 = M G 2 + NG 3 ; (14.4) 



and inverse one given by 



G 2 = N F 2 - (3 F 3 , 
G 3 = -M F 2 + aF 3 . (14.5) 



Combining equations in (|14.3p . we get 



„ d „ , T „ T . J 1 m 2 . _ 

2ti , G 2 -N F 2 + N (-— — + F 3 - 

dy y(i + y) i + y y(i + y) 

fcm „ „ . u 2 k 2 , „ 

y(i + y) y y(i + y) 

2 ^ G + M ^ m F M ( ^ + ^ ) F + 

dy y(i + y) i + y y(i + y) 

km _ . w 2 A; 2 . _ 

y(i + y) y y(i + y) 



(14.6) 
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Expressing F 2 ,F 3 through G 2 ,G 3 according to (|14.4|) : 



2lo— h 

dy 



-(Na + pM) 
-2Np 



km 



y(i + y) 

km 



+ NM 



-to 2 y + m 2 

y(i + y) 



Pa 



u 2 {l + y)-k 2 

y(i + y) 



G 2 + 



+ N- 



2 —ojy + m 



\i + y)-k 2 
y(i - y) 



G 3 = 0, (14.7) 



2uj— h 

ay 



{MP + aN) 



km 



+ 



2Ma 



y(i + y) 

km 



NM-" 2y + m2 +Pa ^ + v)-^' 



y(i + y) 



M 



., — u; 2 y + m 2 9 
2 y +a 2 



y(i + y) y(i + y) 

Let us detail the coefficients at G 3 and G 2 : 



y(i + y) 

w 2 (l +y) - A: 2 n 

y(i + y) 



G 3 + 



G 2 = , (14.8) 



-2A7? 



/cm 



+ iV' 



2 — u> 2 y + m 2 

y(i + y) 



/3 ; 



2 {l + y)-k 2 



y(i - y) 



G, 



y(i + y) 

' [ -2£;m Nfi - uj 2 {N 2 + p 2 ) y + N 2 m 2 + p 2 k 2 } , 



2Ma 



y(i + y) 

km 



y(i + y) 
l 



, , 9 — w 2 y + m 2 9 w 2 (l + y) — k 2 
M ^ — + a 2 v 



y(i + y) 



y(i + y) 



Go 



[ 2km Ma + u?(M 2 + a 1 ) y - M' A m 2 - a 2 k 



2 _2 2 ,2 



y(i + y) 

Redundant singularities will be excluded if 



N 2 + P 2 = , M 2 + a 2 = , 



With additional assumption of unitarity for transformations (|14.4p and (|14.5p 

a G 2 + P G 3 = cos A G 2 + i sin A G 3 , 
M G 2 + N G 3 = i sin A G 2 + cos A G 3 , 



we get 



Na + PM = cos 2A , 2NP = isin2A, 2Ma = isin2A 



aP = i sin A cos A = — sin 2 A , 

N 2 = cos 2 A , /3 2 = - sin 2 A 
correspondingly eqs. (|14.7|) and (|14.8|) give 

km 



9 d ° 2 a 
2cj— h 

dy 

-i sin 2A 



cos 2 A 



y(i + y) 2 



+ - sin 2 A 



NM = i sin AcosA = - sin 2^4 , 

2 

M 2 = - sin 2 A , a 2 = cos 2 A ; 
-a; 2 ?/ + m 2 — uj 2 (1 + y) + k 2 



km 

y(i + y) 



y(i + y) 

,2 c 



2 -u^y + m 2 . 2 w 2 (l + y)-A; 2 
+ cos A ; h sm A 



y(i + y) 
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y(i + y) 



G 2 + 



G, = 



(14.9) 



(14.10) 



2ui— h 

ay 



cos 2 A 



km 



y(i + y) 2 



- sin 2 A 



-u) 2 y + m 2 — oj 2 (1 + y) + k 2 



y(i + y) 



G 3 + 



+ 



i sin 2,4 — r + sin A — + cos 2 A 



2/(1 + y) vO- + y) 

With additional requirement 

1 



+ y) 



G 2 = o 



A = vr/4, cos 2 vl = sin 2 yl = - , sin2A = l, cos2,4 = 0; (14.11) 



the system becomes much more simple 
dG 2 .-uj 2 y + m 2 - uj 2 {l + y) + k 2 



2uj—^ +i 
dy 



2y(l + y) 



-2ikm - u) 2 y + m 2 + uj 2 (\ + y) - k 2 

G2 + mrv) ° 3 = 



dG% .—uj 2 y + m 2 —uj 2 (\ + y) + k 2 2ikm — uo 2 y + m 2 + uj 2 (1 + y) — k 2 

2ul — ; I — — ; G3 H - G 2 = 



dy 



2y(l + y) 



2y(l + y) 



or 



d .ui 2 y + ui 2 {l + y)-m 2 -k 2 \ ui 2 + (m - ik) 2 

2iO— l ; : Go H ; : Cxa = U , 

dy 2y(l + y) ^ 2y(l + y) 



n d ,uj 2 y + uj 2 (l+y)-m 2 -k 2 \ lu 2 + (m + ik) 2 

2UJ- h l ; : H ; r Go = . 

dy 2y(l + y) J 3 2y(l+y) 

Prom whence it follows that 

2y(l + y) riG 2 . Vy + u; 2 (l + y) - m 2 - k 2 

G3 — — ZU ; / TTo - ; r » 9 ; ; rj^ G2, 



(14.12) 



d 2 G 



uj 2 + (m — ifc) 2 w 2 + (m — ik)' 



2 / 9 m 2 m 2 + fc 2 



%(i + ^ + 4 ( i + 2 S ,^ + ^ + ^- !/(i+!/) t 



After changing the variable y to — y eq. (|14,13h reads 

d 2 G 2 „ ,dGi 



G 2 = . 



(14.13) 



4y(l-y)^-/ + 4(l-2y) 



dy 



-2iu + uj 2 + 



m 2 m 2 + k 2 



Making substitution 



y(i-y) l-y 

G 2 = y A (l-y) B G(y), 



Go = 



(14.14) 
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we get 

d 2 G dG 
4y(l - y)-^r + 4[1 + 2A - (2A + 2B + 1 + l)y] — 

AA 2 -m 2 4B 2 + k 2 ' 



-w(2» - w) + 4M + 5)M + B + 1) 

y 1 - y 

With additional requirements 

m = +2A , k = ±2iB . 

it becomes an equation of hypergeometric type 

fin An 

4y(l - v)~rj + 4[1 + 2A — (2A + 2B + 1 + l)y] — 
dy dy 

- [-u(2i - to) + A(A + B)(A + B + 1)] G = 

7=1 + 2A , a + /? = 2,4 + 25 + 1, 
-w(2i - w) + 4(A + + B + 1) 



G 



that is 



a/3 



a = A + B-y, /3 = A + B + 1 + — . 

In the same manner, from eqs. (|14. 12[) it follows 

2y(l + y) dG 3 . u; 2 y + lo 2 (1 + y) - m 2 - k 2 
^2 — ~~ liO— -. tttk — ; 1 n 7 Tv5 "3 

or + (m + ik) dy <jj A + (m + 7/cj z 

d / 2y(l+y) rfG 3 . ^ 2 y + ^ 2 (1 + y) - ^ 2 - k 2 

lid— I — liO—^ TTTK—, 1 n / TT9 "3 

dy \ uj* + (m + iky dy <jj £ + (m + ik) A 

i m 2 y + LU 2 (l+y)-m 2 -k 2 ^ 2y(l + y) dG 3 _ 
2y(l + y) a; 2 + (m + z/c) 2 

,uj 2 y + uj 2 (\ + y) - m 2 - k 2 to 2 + (m - ik) 2 

w 2 + (m + ifc) 2 3) + 2y(l + y) 3 



or 



%(1+9) ^ +4(1+2 ^ + 



/ 9 m 2 m 2 + fc 2 \ „ 



After changing the variable y to —y 



4*1-,)** +4(1-*)^ 



77) 



2 



to 2 + /c 2 



2io; + w 2 + - - ' - ) G 3 = 

y(i-y) i-y 
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and with the substitution 

G 3 = y A (l-y) B G(y), 

we get 



d 2 G rIC 
%(1 - y)-^ + 4[1 + 2A - (2A + 2B + 1 + l)y] — 



n . , „ . AA 2 -m 2 \B 2 + k 2 
u{2i + uj) + A(A + B)(A + B + 1) ^— 

With the help of additional restriction 

m = ±2A , k = ±2iB . 
the latter reads as an equation of hyper geometric type 



G = 



±± + 4 [1 + 2A - (2A + 2B + 1 + l)y] ^ - 

«y «y 

- [u(2i + w) + A(A + B)(A + 5 + 1)] G = , 



c = 1 + 2,4 , a + 6 = 2A + 2B + 1 , 
u(2i + Lo)+4(A + B)(A + B + 1) 



that is 



ab 



IU1 IU1 

a = A + B + l- — , b = A + B + -. 



Let us find a relative factor in two functions G2 and G 3 . Starting from 

G 2 = M 2 yH/2 n _ y \|fc|/2i F (_ n; n + 1+ I m I I m I +1; y) 

i 



M 2 y^' 2 (l - y)( a+b - c V 2 F(a,b,c,y) 



Then 



G 3 = M 3 y H/2 * (1 - y) W2i F{-n + 1, n+ | m | | m | +1; y) = 

z 

= M 3 yM^l - y)^ a+b - c ^ 2 F{a + 1, 6 - 1, c, y) . 

and allowing for 

G 3 [(m - ifc) 2 + a; 2 ] = -4a; y(l -y)^ + i[-m 2 - k 2 + w 2 (l - 2y)] G 2 

dy 
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we arrive at 

M 3 

(m — ik — ico) (m — ik + ioj) -rj~ Fziv) = 

= -4w [ ^ (1 - y) F 2 (y) -^y F 2 (y) + 
2 2i 

+y(l -y) — F 2 (y) ] + i[-m 2 - k 2 + uo 2 {l - 2y)] F 2 (y) . (14.29) 
ay 

To find the relative factor it sufficient to consider the latter at the point y = which results 

in 

(m — ik — iuS) (m — ik + iu) — - = —2uj I m I —in? — ik 2 + iu; 2 . 
v ' y ' M 2 ' 



or 



that is 



M 3 

i(—iuj + m — ik) (iuj + m — ik) —— = (— iuj+ \ m \ —ik) \ m \ +ik) 

M 2 



M 2 = iM {—iuo + m — ik)(iu + m — ik) , 
M 3 = M (-iuj+ J m I -ik)(iuj+ \ m \ +ik) . (14.30) 

These relations become more simple when separating regions for m: 

m > , M 2 = iM(iuj - ik + m) , M 3 = M(iuj + ik + m) ; 

m < , M 2 = iM(iuj - ik + m) , M 3 = M(w + ik + m) ; 

m = , M 2 = M(k - u) , M 3 = iM(k + w) . 

15 Discussion: on relation to other formalisms 
in Maxwell theory 

There exist close relation between the above used covariant technique in the complex form 
of Riemann-Silberstein-Majorana-Oppenheimer and spinor form of Maxwell theory in general 
relativity [18] mainly used in the form of Newman- Penrose formalism of isotropic tetrad 



ia a (x) [ d/dx a + E a (x) ® I + I ® E Q (z) ] £(x) = - j(x) , 
ia a (x) [ d/dx a + E a (x) ® I + J® t a (x) } rj(x) = - j(x) , 

^F af3 (x) = - j a (x) , e Q ^V pF pa (x) = ; (15.1) 

ia a (x) [ d/dx a + E a (x) 01 + 1® S a (x) ] H(x) = £ , 

ia a (x) [ d/dx a + T, a (x) 01 + 1® S a (s) ] A(x) = rj , 

V a A a = , V a A p - VpA a = F a(] ; (15.2) 
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where spinor connections by Infeld - van der Vaerden are used 




), a a (x)=a a ef a) (x), 
» = \ S ab e[ x) V„(e (6)/3 ) , 

S~ a 1 ( a -b _b -a \ . 

= — { <T (J — (J (J ) ; 



(15.3) 



and electromagnetic bi-spinor and electric source are given by 



H r/ 



+ 



S mn a 2 E, 
-i a 1 a 2 Ai 



mn 



+i a 1 a 1 Ai 
- t mn a 2 F n 



mn 



j(x) = i <j k a 2 j k (x) 



j(x) = -i a k a 2 j k (x) . 



(15.4) 



Used in the paper complex 3-vector approach seems to be a good alternative to the other 
possible techniques. Let us summarize the content of the paper again. 

Complex formalism of Riemann - Silberstein - Majorana - Oppenheimer in Maxwell electro- 
dynamics is extended to the case of arbitrary pseudo-Riemannian space - time in accordance with 
the tetrad recipe of Tetrode - Weyl - Fock - Ivanenko. In this approach, the Maxwell equations 
are solved exactly on the background of static cosmological Einstein model, parameterized by 
special cylindrical coordinates and realized as a Riemann space of constant positive curvature. 
A discrete frequency spectrum for electromagnetic modes depending on the curvature radius of 
space and three parameters is found, and corresponding basis electromagnetic solutions have 
been constructed explicitly. In the case of elliptical model a part of the constructed solutions 
should be rejected by continuity considerations. 

Similar treatment is given for Maxwell equations in hyperbolic Lobachevsky model, the 
complete basis of electromagnetic solutions in corresponding cylindrical coordinates has been 
constructed as well, no quantization of frequencies of electromagnetic modes arises. 

16 Supplement: on the use of matrix complex form of the 
Maxwell equations construct electromagnetic solutions from 
scalar ones in Minkowski space-time 

The above matrix form of Maxwell theory : 



being applied to the flat Minkowski model, permits us to develop a simple method of finding 
solutions of Maxwell equations on the base of known solutions of the scalar massless equation 
by Klein - Fock - Gordon. Indeed, in virtue of the above commutative relations we have an 
operator identity 







(-ido + a J dj)V = , 



* = 



(16.1) 



(—ido + a l d\ + a 2 d2 + oi 



3 d 3 ) {-id ~ a 1 8 l - a 2 d 2 - a 3 8 3 ) = + d\ + d\ + d\) . 
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Therefore, taking any special scalar solution one can immediately construct four solutions of the 
Maxwell equation: 



(ido + a l di + a 2 d 2 + a 3 d 3 ) $(x) = 

= {I' ,* 1 ,* 2 ,* 3 } 



id $ 




d 2 <S> 


a 3 $ 


-di<s> 






<9 2 $ 




d 3 $ 






-d 3 $ 


-d 2 $ 


9i$ 





(16.2) 



Thus, we have four formal solutions of the free Maxwell equations (let F a (x) = d a &(x)): 



iF 


Fi 


F 2 


F 3 


-Fi 


iF 


-F 3 


F 2 


-F 2 


F 3 


iF 


-Fx 


-F 3 


-F 2 


F x 


iF 



(16.3) 



In general, the function 3>(x) is a complex- valued. Relationship defining all possible solutions 
has the structure 



A *° + Ai^ 1 + A 2 * 2 + A 3 * 2 





B + icB 



116-4) 



Substituting these expression into Maxwell equations and performing calculations needed we 
can derive four relations (for more details see |23j): 



[ - A d + % Xj dj]F c = 0. 



16.5) 



These equations should be analyzed instead of Maxwell equations to construct solutions of them 
in any coordinate system. 

It seems reasonable to expect further developments in this matrix based approach to Maxwell 
theory, as a possible base to explore general method to separate the variables for Maxwell 
equations in different coordinates (see 
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